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We calculate the effect of variation in the light-current quark mass, 



on standard big bang 



nucleosynthesis. A change in m q at during the era of nucleosynthesis affects nuclear reaction rates, 
and hence primordial abundances, via changes the binding energies of light nuclei. It is found 
that a relative variation of 8m q /m q = 0.016 ± 0.005 provides better agreement between observed 
primordial abundances and those predicted by theory. This is largely due to resolution of the existing 
discrepancies for 7 Li. However this method ignores possible changes in the position of resonances in 
nuclear reactions. The predicted 7 Li abundance has a strong dependence on the cross-section of the 
resonant reactions 3 He (d,p) 4 He and t (d, n) 4 He. We show that changes in m q at the time of BBN 
could shift the position of these resonances away from the Gamow window and lead to an increased 
production of 7 Li, exacerbating the lithium problem. 
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I. INTRODUCTION 

Measurements of the primordial baryon-to-photon ra- 
tio r\ from the cosmic microwave background from 
WMAP [l[, coupled with precise measurements of the 
neutron half-life [2| , have made big bang nucleosynthesis 
(BBN) an essentially parameter- free theory 0, [1, H| . In 
this paradigm excellent agreement has been obtained be- 
tween predicted and observed abundances of deuterium 
and 4 He (see, e.g. the Particle Data Group review [2| and 
references therin). However there is some disagreement 
for 7 Li, the only other element for which the abundance 
has been measured to an accuracy at which fruitful com- 
parison with theory can be made. While the "lithium 
problem" has been known for some time, it has been ex- 
acerbated by recent measurements of the 3 He(a, 7) 7 Be 
reaction Standard BBN theory with r\ provided 

by WMAP 5 overproduces 7 Li by a factor of 2.4 - 4.3 
(around 4 - ha) 

One possible solution to the lithium problem is that 
the physical constants of the early Universe may have 
been slightly different. In fact, such variations in the 
physical laws can be well-motivated theoretically in an 
expanding Universe; see @] for a review. Ref. 0] consid- 
ered variation of the deuterium binding energy B d during 
primordial nucleosynthesis. BBN has a high sensitivity 
to B d since its value determines the temperature at which 
deuterium can withstand photo-disintegration and hence 
the time at which nucleosynthesis begins. Their best- 
fit result AB d /B d = -0.019 ±0.005 resolved then-extant 
discrepancies between theory and observation in both Li 
and rj (or alternatively, in 7 Li and 4 He with 77 fixed by 
WMAP). 

More recently, Ref. Q examined the response of BBN 
to variation of several physical parameters, including 
binding energies, in a linear approximation. These were 
coupled with calculated dependences of binding energies 
on m q in 9], which found that the 7 Li abundance dis- 
crepancy could be resolved by a variation in light-quark 



mass of Sm q /m q = 0.013 ± 0.002. Crucially, the 4 He and 
d abundances were found to be relatively insensitive to 
m q and so the existing agreement between theory and 
observation in these elements was maintained. 

In this paper we re-examine the dependence of light- 
element production on variation of the dimensionless pa- 
rameter X q — m q / Aqcd where m q is the light-quark 
mass and Aqcd is the pole in the running strong- 
coupling constant. We follow @ and assume that Aqcd 
is constant, calculating the dependence on the small pa- 
rameter m q . This is not an approximation. Rather it 
only means that we measure all dimensions (m q , cross 
sections, etc) in units of Aqcd- Therefore 5m q /m q 
should be understood as SX q /X q . We take into account 
several effects that were not previously considered, most 
importantly the nonlinear dependence on m q and varia- 
tion of resonance positions. 



II. VARIATION OF BINDING ENERGIES 

The energy released in each reaction, Q, is determined 
by the masses of the reactants and products, which in 
turn are determined by the nuclear binding energies. As 
noted in [8| , the Q- values affect the forward (exothermic) 
reaction rates via phase space and radiative emission fac- 
tors. For radiative capture reactions at low energy E the 
Q-dependence is 



a(E) cx E* 



{Q + Ef 



(1) 



For low-energy reactions with two nucleons in the exit 
channel the dependence is proportional to the outgoing 
channel velocity, v ~ (Q+E) 1 / 2 . When the outgoing par- 
ticles are charged, the Gamow factor of the exit channel 
can also contribute: 



a(E) ~(Q + E) 



1/2 -y/B e /(Q+E) 



(2) 



The Gamow factor appears because of the Coulomb bar- 
rier to the reaction; E g — 2ir 2 ZfZ^a 2 fie 2 where a is the 



2 



TABLE I: 9 In Ya/dlnBo, the dependence of nuclear abun- 
dances, Y a , on deuterium binding energy under different 
assumptions: 

1. Variation of virtual level not considered, (av) ~ Q 5 ^ 2 . Q 
changed only for p (n, 7) d. 

2. Variation of virtual level not considered- effect of Bd 
included in all reactions (similar to theory of |8|). 

3. p(n,"/)d changed according to ([5]), including variation 
of the virtual level; effect of Bd on other reactions ignored 
(similar to theory of Q). 

4. p (n, 7) d changed according to ((5]); effect of Bd included in 
all reactions. 



Method 


d 


J He 


4 He 


b Li 


; Li 


f. 


-4.04 


-1.75 


0.68 


-3.17 


10.59 


2. 


-2.91 


-2.08 


0.67 


-6.58 


9.41 


3. 


-5.f2 


-1.29 


0.70 


-4.23 


17.99 


4. 


-4.00 


-1.62 


0.69 


-7.64 


16.81 



fine-structure constant, Z\ and Z2 are the charge num- 
bers of the products, and fi is the reduced mass of the 
products. At BBN temperatures we can usually assume 
that E <C Q- Expanding in Q, 




a = (T 



and we see that the Gamow term in @ is generally 
small (it was neglected in [8j). However it can be im- 
portant for some reactions, for example in 7 Be (n,p) 7 Li, 
\J Eg/Q — 2.17, i.e. it triples the effect of SQ on the 
reaction rate. 

The reverse reaction rates are simply related to the for- 
ward rates via statistical factors. From detailed balance 
one finds 



'rev 



-Q/t 



{(TV 



fwd 



(4) 



and we see that the reverse reactions also provide sensi- 
tivity to Q. 

An exception to the rule (jTJ) is found in the reaction 
p (n, 7) d, an important reaction because d is a precursor 
to all further nucleosynthesis. This reaction is sensitive 
not only to Q but also to the position of the virtual level 
with energy e v — 0.07 MeV. The sensitivity of this reac- 
tion to Q was calculated in Q 



{(TV 




(5) 



Note that 8, 9] did not take variation of the virtual level 
into account. In Table[T]we show the linear dependence of 
abundances on the deuterium binding energy with differ- 
ent theories of variation. It shows the effect of variation 
of the virtual level, as well as the effect of including Bd 
variation on other Q-values and reaction rates. 



We denote the sensitivity of nuclear binding energies 
to the light-current quark mass m q by 



K 



SE/E 
Sm q /m q 



(G) 



Values of K for several light nuclei were presented in 
Refs. We use the "best values" from these pa- 

pers, given by the AV18+UIX nuclear Hamiltonians, 
with hadron mass variations calculated in terms of the 
m n using the Dyson-Schwinger equation calculation of 
[111 ] . From these one calculates the m g -dependence of the 
Q values, and therefore the reaction rates, and therefore 
the primordial abundances of light elements in BBN. 

In Fig. Q] we present our predicted values of 4 He, 
d, and 7 Li with different values of light-quark mass. 
Details of the calculations and explanation of obser- 
vational abundances are presented in the appendices. 
Comparing the observed and predicted abundances from 
the figures we obtain for 4 He, d, and 7 Li respec- 
tively, Sm q /m q = -0.002 ± 0.037, 0.012 ±0.011, and 
0.018 ± 0.006. The three data sets are therefore consis- 
tent, with weighted mean 



5m q /m q = 0.016 ± 0.005 



(7) 



(3) It is seen that the 4 He abundance has a low sensitivity 
to m q ; furthermore we show in Fig. [1] the conservative 
observational error bounds provided by [12fl . Therefore, 
it is worth pointing out that the more tightly constrained 
abundance, Y p = 0.2477 ± 0.0029 [Hj], is also consistent 
with the variation ([7]). It is clear from Fig. [T]that tak- 
ing into account the nonlinear dependence of BBN abun- 
dances on m q is important, particularly for 7 Li. In fact, 
if we assume a linear response, as was done in [1,01, we 
instead obtain 5rria/m q = 0.014 ± 0.002. 

As noted in 9], to take into account uncertain- 
ties in the theoretically derived quantities K (Equa- 
tion [6|) the final result ([7|) should be interpreted as 
5m q /m q — k ■ (0.016 ± 0.005) where k ~ 1 and the ac- 
curacy in k is approximately a factor of two. 



III. RESONANCES 

Of the most important reactions in BBN, the mirror 
reactions 



3 He (d,p) 4 He 
t(d,n) 4 Ue 



(Reaction 1) 
(Reaction 2) 



are the only reactions where the cross-section is dom- 
inated by a fairly narrow resonance. Therefore, one 
can hope for sensitivity of primordial abundances to the 
position of these resonances. (Note that the reaction 
7 Be (n,p) 7 Li is also dominated by a near-threshold res- 
onance, however in this case the resonance is a rather 
broad and hence strong sensitivity can hardly be ex- 
pected.) 
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FIG. 1: Calculated 4 He, d, and 7 Li abundances vs. rela- 
tive change in light quark mass m q / Aqcd (solid lines). The 
ranges showed by the dashed lines are la errors in the theory, 
assuming the relative errors are constant (i.e. these do not 
take into account any error in the K factors of Eq. [6}. The 
shaded areas show la ranges of observed abundances (details 
in Appendix |A)| , 



Both of these reactions have the cross-sections with the 
general form 



a(E) = 



P(E) 



E {E-E r ) 2 +Tf/4 



(8) 



where E g is the Gamow energy of the reactants, E r 
and F r are resonance parameters, and P(E) is a poly- 
nomial chosen to fit the measured reaction cross-section. 
In this work we use the cross-section fits of Ref. [l4j |. 



which 

,(2) 



(1) 



0.183 MeV, n ' = 0.256 MeV and 
E?> = 0.0482 MeV, T^ 2) = 0.0806 MeV for reactions 1 
and 2, respectively. 

Consider modification of the resonance positions, 
E r — > E r + SE r , due to a variation of the fundamental 
constant m q . Reaction 1 will be affected in the following 
way. The resonance is an excited state of 5 Li; that is, a 
compound nucleus with three protons and two neutrons: 
we call this state 5 Li . Similarly there is a state 5 He for 
reaction 2. Then 



EM 



Es 



Li* " 
He* 



E. 



^He 



E d 



E± — Ed 



(9) 
(10) 



and so E, Lr = -9.76 MeV and £ sHe * = -10.66 MeV. 
The change in the resonance position due to a variation 



in m q is therefore 



5E<? 



6E 5U * - 5E 3Re - SE d 



(11) 



{K*U' E *W - ^He^He - K d E 4 ) — £2) 



with the K defined by ©. 

Changes to the cross-section of reaction 1 affects the 
primordial abundances of 3 He and 7 Be, while changes 
in reaction 2 affect abundances of t and 7 Li. Since t and 
3 He are not well constrained observationally, we choose to 
focus on Li. In Fig. [2] we present Li abundance against 
variation of light quark mass Sm q /m q at rj = 6.23 x 10 -10 , 
the WMAP5 value. For such a value of rj, the majority 
of 7 Li is created as 7 Bc (which /3-captures to 7 Li) via the 
reaction 3 He ( 4 He, 7) 7 Be. 

We need to find Kr,^* (and similarly if B jj e *). One as- 
sumption is that the mass-energy of the resonance varies 
with the mass-energy in the incoming channel Q ; in this 
case the resonance does not shift. This assumption cor- 
responds to -ft^Li* = —1.54 and ^ 5 jje* = ~ 1-44. It 
corresponds to the solid line in Fig. O 

A more reasonable guess is to assume that the vari- 
ation of the resonant state 5 Li will be approximately 
the same as that of the ground state 5 Li. This can be 
seen by considering the resonance and the ground state 
configurations as residing in the same potential. The sen- 
sitivity of the ground state 5 He to m q has been calculated 
^ 5 He = —1-24 [nj; K 5 ^ was not calculated explicitly, 
but its value will be very close to that of 5 He. Our as- 
sumption of equal variation of the ground and excited 
state then gives 



K 5U * = -3.35 



K 



J>He* 



3.19 



(13) 
(14) 



This assumption corresponds to the dashed line in Fig. [5] 
The equal-variation assumption in the previous para- 
graph represents an upper limit on the relationship be- 
tween the ground and excited state. In reality the 
potential-dependence of the states may be different, in 
which case the shift of the 5 Li (or 5 He) resonance may be 
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FIG. 2: Calculated 7 Li abundance vs. relative change in light 
quark mass m q / Aqcd ■ Solid line: no shifts in resonance po- 
sitions included (same as solid line in Fig. [TJ; dashed line: 
resonance shifts according to assumption that resonant state 
varies as much as the ground state fequations ll3l andll4 [l : dot- 
dashed line: an averaged value of resonance-position sensitiv- 
ity used. The shaded area shows the la range of abundances. 



smaller than the shift of the ground state. On the other 
hand a minimum value of K for the resonance states 
is that of the ground state, Ks\jy* = K 5 = —1.24. 
A reasonable, conservative, estimate is to take the av- 
erage of these extremal values: K 5 ^* = —2.29 and 
^ 5 He* = —2.21; this is the dot-dashed line in Fig. [5] 
Ultimately however, we require a nuclear calculation of 
sensitivity, of the kind presented in Refs. (9l, [Io|. 

The effect of SEP on BBN can be understood in 
the following way. When the cross-section is convolved 
with a Maxwellian distribution, the exponential term 
gives rise to the "Gamow window" at energy Eo/E g — 
{kT/2E g ) 2 / 3 . This reaction is most active at kT w 0.07 
MeV, at which time the Gamow window is at Eq = 
0.180 MeV. This is remarkably close to the resonance 
energy for this reaction E r = 0.183 MeV. Therefore 
movement of the resonance position in either direction 
will reduce the cross-section for this reaction at the rel- 
evant temperatures. In turn this reduces the amount of 
3 He that is destroyed via reaction 1, leaving more to react 
with 4 He to produce 7 Be. On the other hand the effect 

of this reaction on d and 4 He abundances is minimal. 

(2) 

The effect of 5E). is very similar: it reduces the 
amount of t destroyed in reaction 2, leaving more tri- 
tium to react with 4 He to produce 7 Li directly. De- 
spite this production channel being suppressed at high 
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FIG. 3: Calculated 3 He abundance vs. relative change in 
light quark mass m q / Kqcd- Solid line: no shifts in resonance 
positions included; dashed line: resonance shifts according to 
assumption that resonant state varies as much as the ground 
state (equations [13] and I14[) : dot-dashed line: an averaged 
value of resonance-position sensitivity used. 

1], the effect of SE^ is still important for 7 Li produc- 
tion because the relative effect of the variation is larger: 
5E^ ] /T ( r ] > 5eP/T^\ The trends seen in Fig. [2] are 
the same even at low 77 since both reaction pathways be- 
have in much the same way to variation in m q . 

From Fig. [5] we see that taking shifts in the resonance 
positions into account can destroy the agreement between 
theory and observation previously obtained by varying 
Tdq . In the case where the shifts in the ground and reso- 
nant states vary by the same amount (dashed line), the 
7 Li discrepancy actually gets worse with variation in light 
quark mass. On the other hand the milder "averaged K" 
response (dot-dashed line) still significantly challenges 
the conclusions of Section [TTJ It is not appropriate to 
directly compare primordial 3 He abundances with obser- 
vations because of the complexity of the stellar evolution 
of this isotope [l5| , however we note that primordial 3 He 
production could also be greatly increased by movement 
of these resonances (Fig. [3]). 



IV. CONCLUSION 

We have shown in Section |TT] that a variation in the 
light-quark mass during the era of big-bang nucleosyn- 
thesis of Sm q /m q — 0.016 ± 0.005 provides better agree- 
ment between theory and the observed primordial abun- 
dances. This is largely because it resolves the existing 
disagreement in 7 Li abundances 0]. 

However this conclusion is threatened when movement 
of the resonance positions in the reactions 3 He (d, p) 4 He 
and t [d, n) 4 He is taken into account. These reactions 
strongly affect 7 Li production during BBN; furthermore 
they are already "on resonance" meaning that movement 
of the resonance position in either direction increases Li. 
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Our estimates suggest that the 5 He* and 5 Li* resonances 
may be very sensitive to variation of m q / Aqcd- There- 
fore it is very important that the sensitivity of these res- 
onances to fundamental constants be studied in more de- 
tail using nuclear models. 

This work is supported by the Australian Research 
Council. JCB thanks Daniel Wolf Savin for hospitality 
during the early stages of this work. 



APPENDIX A: OBSERVATIONAL 
ABUNDANCES 

Our observational abundances largely follow the rec- 
ommendations of the Particle Data Group review 
The deuterium abundances are derived from several stud- 
ies of isotope-shifted Ly-a spectra in quasar absorption 
systems. Combined these give 

d/H = (2.84 ±0.26) x 1CT 5 

where the errors have been increased to account for the 
scatter between different systems. 

4 He is observed in H II regions of low-metallicity dwarf 
galaxies. A very conservative estimate of observed He 
abundance comes from [l^] 

Y p = 0.249 ±0.009 . 

The error here is significantly larger than other extrapo- 
lations to zero metallicity, e.g. [l3,[l(|. 

Primordial Li abundance is determined from metal- 
poor Pop II stars in our galaxy. Lithium abundance 
does not vary over many orders of magnitude of metal- 
licity in such stars; this is the Spite plateau (T3l - Recent 



studies give abundances of (1.1 — 1.2 ± 0.1) x 10 10 
(1.26 ± 0.26) x 10~ 10 [19], and (1.1 - 1.5) x 10~ 10 
Significantly higher results were obtained with different 
methods of obtaining effective temperature of the stars 
[2~fl ] since the derived lithium abundance is very sensitive 



to temperature. However no evidence for high tempera- 
tures was found in the studies [HI, [2(| ■ 

On the other hand measurements of 7 Li abundance 
in the globular cluster NGC 6397 give values of (2.19 ± 
0.28) x 10" 10 [2|L (1.91 ±0.44) x 10~ 10 [H|, and (1.69± 
0.27) x 10" 10 ' 24]. The M 92 globular cluster yields a 
value of (2.29 ± 0.94) x lO" 10 [2j|. For a more detailed 
review and discussion of systematics, see, e.g. (2L 14. Il8j. 

In this paper we use the conservative range of 



7 Li/H = (1.5 ±0.5) x 10 



-10 



which was also adopted in [8j, although we note that 
some of the studies listed above give ranges as high as 
2.5 x 10- 10 . 



APPENDIX B: COMPUTER CODE, REACTION 
RATES, AND THEORETICAL UNCERTAINTIES 

In this work we calculate BBN abundances using a 
modified Kawano code [26| , with updated reactions from 
the NACRE collaboration [13] • For the most important 
reactions we use the cross-section fits in [3] ■ The excep- 
tions are p(n,^)d where we use the calculations of [28j . 
and the recently measured reaction 3 Hc ( 4 He, 7) 7 Be for 
which we use the fits provided in [5j. We have not cal- 
culated errors in the theoretical prediction; instead we 
simply take the relative errors from the recent calcula- 
tions of Cyburt, Fields and Olive [4[ which use a very 
similar reaction network and the latest physical data: 



d/H 
7 Li/H 



0.2486 ±0.0002 
(2.49 ±0.17) x 10" 5 
(5.24±°;£) x 10- 



-10 



(Bl) 
(B2) 
(B3) 



These values compare well with the results from 
our code of Y p = 0.2486, d/H = 2.53 x 10~ 5 , and 
7 Li/H = 5.05 x 10~ 10 . 
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